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All stationary solutions to the one-dimensional nonlinear 
Schrodinger equation under box or periodic boundary condi- 
tions are presented in analytic form for the case of attractive 
nonlinearity. A companion paper has treated the repulsive 
case. Our solutions take the form of bounded, quantized, 
stationary trains of bright solitons. Among them are two 
uniquely nonlinear classes of nodeless solutions, whose prop- 
erties and physical meaning are discussed in detail. The full 
set of symmetry-breaking stationary states are described by 
the C n character tables from the theory of point groups. We 
make experimental predictions for the Bose-Einstein conden- 
sate and show that, though these are the analog of some of the 
simplest problems in linear quantum mechanics, nonlinearity 
introduces new and surprising phenomena. 



I. INTRODUCTION 

In a recent experiment the dilute-gas, attractive, Bose- 
Einstein condensate (BEC) was created for the first 
time for lithmrrJiTH. As predicted by the nonlinear 
SchrodingeiLm equation (NLSE) for three dimensionsQ, 
the condensate collapsed when the number of particles 
became large. However, in one or quasi-one dimension 
no collapse is predicted!!!, and in one dimension the NLSE 
has application in fiber opticso as well as other 

fieldsffM 

In this article we present stationary solutions to the 
one-dimensional NLSE for attractive nonlinearity under 
periodic and box boundary conditions. In the proceed- 
ing articles we solved the analogous problem for repul- 
sive nonlinearity. Although we will place our results in 
the context of attractive atomic interactions in the BEC, 
they are equally applicable to wave phenomena, in many 
physical situations, as for example ring laserstfl. 

We define the BEC to be in the quasi-one-dimensional 
regime when its transverse dimensions are on the order 
of its healing length and its longitudinal dimension is 
much longer than its transverse ones. In this case the 
ID limit of the 3D NLSE_is appropriate, rather than a 
true ID mean field theoryHl, as would be the case for a 
transverse dimension on the order of the atomic interac- 
tion length. Under these criteria the condensate is well 



out of the Thomas-Fermi limit; i.e. the kinetic energy 
in the transverse dimensions is very high. It is this high 
kinetic energy which prevents the condensate from col- 
lapsing. We have numerically illustrated the stability of 
the condensate in quasi-lD elsewhereQ. ,— ■ 

As was shown in the previous articled, in the quasi- 
lD regime the solutions are approximately separable, and 
the dimensionless NLSE may be written: 



-\ 2 d 2 - \f(x)\ 2 +V(x)]f(x) = pf(x) 



(1) 



where f(x) is the dimensionless wavefunction describ- 
ing excitations along L, x S [0,1], V(x) is a dimen- 
sionless potential representing the boundary conditions, 
ft = (2m^ 2 /h 2 )fi, and A = £/L oc N^ 1 / 2 , the ratio of the 
healing length to the box/ring length. We remind the 
reader that £ 2 = l/(87rp | a |), where a is the scattering 
length and p is the mean particle density. 

For comparison with experiment we list the conversion 
factor from the unitless /i_to // in /zK. Using the formula 
from the previous articleEj, for 7 Li with p = 10 14 cm -3 
and a = —1.45 nm, the conversion factor is 0.125. Since 
the unitless chemical potentials we find will be on the 
order of -1 to -10 this gives a sense of the energy scale of 
the solutions, on the order of 0.1 to 1.0 pK. We also note 
that throughout our presentation we will use a reasonable 
test scale of A = 1/25 for illustrative purposes. 

As | f(x) \ 2 is the longitudinal portion of the single 
particle density, we require the normalization condition: 



[ dx\f(x)\ 2 =l 
Jo 



(2) 



be satisfied together with the NLSE ([!]) and such bound- 
ary conditions as we will describe below. These are the 
equations we will solve. 



II. BOX BOUNDARY CONDITIONS 

The key concept underpinning the approach taken in 
this paper is that in a region of constant potential the 
NLSE-is integrable. This was utilized in our companion 
papert 3 !, and virtually all the techniques applied there 
to the case of repulsive nonlinearity apply here as well, 
though they lead to the identification of different partic- 
ular solutions. 



1 



Specifically, it was shown in the previous articled! that 
solutions of the stationary NLSE on the unit interval, 
subject to the box boundary conditions 



/(0) = /(l) = 



(3) 



can be expressed, without loss of generality, by purely 
real functions. Such solutions are given in terms of Ja- 
cobian elliptic functions, whose properties were reviewed 
in the previous article. This result holds independent of 
the sign of the nonlinearity. For the case of repulsive 
nonlinearity, it was established that the sn function pro- 
vides all possible solutions. For the attractive case, both 
the cn and dn functions offer solutions. However, the dn 
function has no real zeros, so the cn function provides all 
solutions which satisfy box boundary conditions. 



A. Solutions and spectra 



The most general form of the solution is: 



f(x) = Acn(kx + S | m) 



(4) 



where < m < 1 is the parameter of Jacobian elliptic 
functions, k and S will be determined by the boundary 
conditions below while A and m will be determined by 
substitution of Eq. (||) into the NLSE and by normaliza- 
tion. 

Since the quarter period of cn is the complete elliptic 
integral of the first kind, K(m), and since cn(0 | m) = 
1, we find that k — 2jK(m) and 5 = —K(m), where 
j € {1,2,3,...}. The physical meaning of j — 1 is the 
number of nodes in the cn function. We will give a more 
general meaning to j below. We then solve Eq. ([[]) by 
substituting Eq. (|j), using Jacobian elliptic identities, 
and setting coefficients of equal powers of cn equal. This 
results in equations for the amplitude squared, A 2 and 
the chemical potential, jl : 



A 1 



2m{2jK{m)) 2 X 2 



fi = -{2jK(m)) 2 X 2 (2m - 1) 



(5) 



(6) 



Substituting Eq. (||) into Eq. (||), utilizing Jacobian 
elliptic identities, and noting that the integral over cn 2 
can be defined in multiples of the quarter period K(m), 
we obtain the normalization condition: 



The plot demonstrates that the solutions are unique. 
Note that the normalization condition requires m to be 
much closer to 1 than in the repulsive case, as may be 
seen by comparing F-ig. [l] to the corresponding figure in 
the previous articletia It follows that the numerics of 
the attractive case are more difficult than those of the 
repulsive one, as may be seen in some of the figures. 

These solutions are in one-to-one correspondence with 
those of the ID particle-in-a-box problem in linear quan- 
tum mechanics. Plots of the wavefunction for the ground 
state and the first three excited states are shown in Fig. |^. 
In Fig. H the chemical potential spectrum of this solu- 
tion type is plotted as a function of A~ 2 , the number of 
healing lengths per box length quantity squared. The 
leftmost portion of the plot corresponds to the particle- 
in-a-box limit and the rightmost portion to the bright 
soliton limit. We now discuss these two limits. 



B. Particle-in-a-box limit 

High chemical potential states in which the peaks over- 
lap become particle-in-a-box type solutions, as can be 
seen in Fig. |2jd. This is both the zero density, linear 
limit and the highly-excited-state limit. Mathematically, 
m — > + and cn — > cos. Physically, jX ^> 1. In this limit 
K(m) -> tt(1/2 + m/8 + 0{m 2 )) and m -► l/OA) 2 , so 
that Eq. (|J) becomes: 

A = i VA 2 (l-^ + 0(m 2 )) 



A = jVA 2 (i- 



l 



2j 2 ^ 2 A 2 +<D( -fX 4 ^ 



(9) 



or putting back in the units: 



j 2 Tr 2 h 2 , 3ra 9xN 



fl - 

2A/L 2 v 



A t j 2 n 



L 2 N 2 



(10) 



which clearly converges to the well-known linear quantum 
mechanics particle-in-a-box chemical potential. Note 
that the first order correction is identical to the repul- 
sive case except that it is negative. One may also obtain 
this result from first order perturbation theory, as has 
been shown in the previous articleEa. 



2(2j) 2 X 2 K(m)(E(m) - (1 - m)K(m)) = 1 



(7) 



where E(m) is the complete elliptic integral of the second 
kind. Eq. (||) then becomes: 

f{x) = V2m(2jK(m))Xcn(K(m)(2jx - 1) | m) (8) 

This leaves the chemical potential (||) and the wave- 
function (||) determined up to the parameter m and the 
scale A. In Fig. [I] a graphical solution of Eq. (0) is shown. 



C. Bright soliton limit 

One may add a peak without disturbing another peak, 
provided that adjacent peaks have opposite phase and 
that the overlap between them is exponentially small in 
the ratio of their separation to their healing length. In 
this limit we ought to recover a series of equally spaced 
sech solutions of alternating phase to within exponen- 
tially small factors, as for example is shown in Fig. ^fa 
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and Fig. ||c. Sech solutions are also called bright soli- 
tons. The sech function is the m — > 1~ limit of the cn 
function in the solution, Eq. (0). 

Bright solitons solve the free NLSE. Thus we should 
find that the wavefunction and chemical potential no 
longer depend on the box length L. In this limit the 
solitons must have a different length scale than A, which 
depended on L. The soliton width is proportional to the 
parameter 



A, 



V ; 



8ttA o 



(11) 



where A t as the transverse area of the box. In the soliton 
literatures, 77 is usually set to 1 by renormalizing the 
wavefunction. 

We now consider the limit A — > 0, which corresponds to 
771 — > 1~. Physically, this means that the peaks become 
highly separated and the interaction between them be- 
comes exponentially small. By using Taylor expansions 
in (1 — m) of the complete elliptic integrals in Eq. (^) we 
find that Eq. (p|) becomes: 



1 



16j 2 A 2 



A* = 

while Eq. (||) becomes: 

/(^) = ^372-^(^(2^-1) I m) 
Putting back in the units we find: 



(12) 



(13) 



Note that Eq. (JlJ) is identical in form to that of the 
Rydberg formula for the energy levels of the Hydrogenic 
atom. Thus to leading order we have found stationary 
state energy levels for whieh the quantum number scales 
as j (harmonic oscillator )tU, j 2 (particle-in-a-box), and 
j (Rydberg formula). This is an impressive illustration 
of the complexity of solution types for stationary states 
of the NLSE. 



III. PERIODIC BOUNDARY CONDITIONS 

There are four solution types for periodic bound- 
ary conditions. There are constant amplitude solutions 
which are plane waves; real, antisymmetric, symmetry- 
breaking solutions, similiar to those found in Sec. [nj real, 
symmetric, symmetry-breaking solutions; and a novel 
class of complex symmetry-breaking solutions. The for- 
mer two are in one-to-one correspondence with particle- 
on-a-ring solutions in linear quantum mechanics; the lat- 
ter two are nodeless and are found only in the presence 
of nonlinearity. As the ring is rotationally invariant, the 
symmetry-breaking solutions will have a high degener- 
acy, in analogy with vortices in two dimensional^. The 
periodic boundary conditions are: 



/(0) = /(I) 



no) = /'(i) 



(17) 



(18) 



1 1 



2777, 1677 2 y 



f(x) _ 1 1 1 

VI - 23/2 J V nU W X " 



S(L) I m) 



(14) 



(15) 



where 8{L) is an offset which depends on the box length. 
But as L —t 00 the offset becomes arbitrary, so that we 
can set it to zero. Note that we put back in the units of 
the wavefunction f{x) which we took out in making the 
separation of variables in the quasi-lD approximation. 

As the reader may verify, for the case in which j = 1 
this is indeed the chemical potential and wavefunction of 
the sech solution to the free ID NLSE: 



(16) 



We have found that this limit suffices to calculate 
chemical potentials for which j < 1/(5 A) to better than 
1%. This estimate assumes an overall scale size of ~ 5£ 
per peak. The L-independent chemical potentials for the 
j = 1, j = 2, and j — 3 solutions shown in Fig. ^| satisfy 
this criterion, for example, as do any ground states for a 
healing length smaller than 1/10. 



A. Constant amplitude solutions 

If we assume that r(x) = constant then we obtain con- 
stant amplitude solutions of the form: 



(19) 



where n € {0, ±1,±2, ...}. The amplitude is constrained 
by normalization to be 1. Substituting Eq. (|l^) into 
Eq. (p]) we find the chemical potential: 



p. = -1 + (A2?to) 2 



(20) 



Unlike for the case of repulsive nonlinearity, where the 
ground state on a ring was the constant solution for 
77 = 0, the ground state of the attractive BEC breaks 
symmetry at some scales, so that the constant amplitude 
solutions are a highl y excited state. This will be discussed 
further in Sec. V A . For 77 ^ each solution is two- fold 
degenerate, as n can be either positive or negative. 

Note that these states could also be termed angular 
momentum eigenstates or quantized jiprtices, as for ex- 
ample in the work of Matthews et al£2 
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B. Real symmetry- breaking solutions 

Of the two Jacobian elliptic functions having differ- 
ent physical form and solving the NLSE for attractive 
nonlinearity, only one could be used for the box, the cn 
function. The dn function satisfied the NLSE but, as 
it never vanished, failed to satisfy box boundary condi- 
tions. But on the ring the dn function poses no difficulty, 
because the boundary conditions are periodic. Thus two 
real solution types are available. Looking at the plot 
of the Jacobian-elliptic functions in the appendix to the 
previous articletLia, we observe that the cn function is anti- 
symmetric with respect to translation by the half period 
2K(m), while the dn function is symmetric under such a 
translation. 



1. Analog to box boundary condition solutions 

As we have exchanged the ring for the box, Eq. (Q) is 
the real solution. One simply changes k from 2jK(m) to 
AjK(m) in order to satisfy Eqs. ([It]) and (|l8|), i.e. from 
multiples of the half period to multiples of the whole 
period. The number of nodes will be 2j rather than j-1, 
where j S {1,2,3, ...}. We temporarily keep 5 set to 0. 
But note that, unlike for box boundary conditions, under 
periodic boundary conditions 6 is arbitrary. 

Then all the results from section || hold with the new 
k, by letting j — > 2j in all equations. The energy and 
wavefunction are determined uniquely by graphical solu- 
tion of Fig. |l|. In Fig. ||d and ||d we show the first two 
states. Both the linear quantum mechanics, particle-on- 
a-ring limit and the bright soliton limits are reproduced. 
In the latter the same kind of non-overlapping criterion 
applies as before. 

We found real solutions by setting (5 = 0. If we instead 
let 8 vary arbitrarily, we obtain the degeneracy inherent 
in these symmetry-breaking solutions. The entropy as- 
sociated with a pair of peaks depends logarithmically on 
the box length L, and, since there are approximately A -1 
possible positions for a peak, the entropy ist2Hl3: 



S ~ fceln( 



1 



10jA' 



(21) 



where the factor of 10 comes from 5 for each of the two 
peaks. This is consistent with the non-overlapping crite- 
rion used in obtaining Eq. (J14|). 



2. Modeless solutions 

To find the dn solutions, we follow the same method 
as outlined in Sec. y. The most general solution is: 



subject to the NLSE the normalization (||), and the 
boundary conditions jl7|) and (|l8|). From these we obtain 
the chemical potential, normalization, and amplitude: 



(X= -{2jK{m)Y\ z {2-m) 



2(2j) 2 X 2 K(m)E{m) = 1 



(23) 



(24) 



f(x) = V2(2jK(m))Xdn(K(m)(2jx - 1) | m) (25) 



where j refers to the number of peaks rather than nodes 
in the dn function, since the dn function is nodeless. Note 
that Eqs. (H)-© are only v alid for m ^ 0. The m = 
case is discussed in Sec. IV A . 

The dn function has a period of 2K(m), not AK(m) 
as did the cn function. Plots of the amplitude for the 
first four energy levels are shown in Fig. ||. The lowest 
state, shown in Fig. ||a, appears identical to the lowest 
state in Sec. [0]. Because m is very close to 1 for the test 
scale of A = 1/25, dn ~ sech and cn ~ sech, so that cn ~ 
dn. Had we chosen a scale at which the peaks came near 
to the boundaries, the two ground state solutions would 
have looked quite different. To make clear the extent to 
which m is singular at this scale, the numerical solution 
to Eq. (§4|) for a single peak is (1 - m) ~ O(10" 66 ). 

The particle-on-a-ring limit for the dn solution is a 
plane wave. As m — > + , dn^ 1. In this case the ampli- 
tude is constrained to be 1 by the normalization and the 
chemical potential is fl — — 1. 

In the bright soliton limit Eq. ( p4| ) may be expanded in 
(1 — m) to yield the same result for the chemical potential 
and amplitude as was found for the cn solutions in a box: 



1 1 



2m 16 j 2 rf 



VI 



-^—-Au{{^-)x ~ 5{L) | m) 
23/2 j 7? yy Ajr]' y ' 1 ' 



(26) 



(27) 



where the same criterion as was used for the cn solutions 
in a box may be applied to the validity of the use of the 
limit here. 

As with the cn solutions, we found the nodeless dn so- 
lutions by setting 5 = 0. If we instead let S vary arbitrar- 
ily, we obtain the degeneracy inherent in these symmetry- 
breaking solutions. The entropy associated with a peak 
depends logarithmically on the box length L, and, since 
there are approximately A^ 1 possible positions for a peak, 
the entropy is: 



S r~j fcB m ( 



1 



nV2j\ 



(28) 



f{x) = Adn(kx + 5 | m) 



(22) where the factor of ny/2 will be explained in Sec. IV A . 
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3. Energy splittings 

Because the dn solutions are symmetric and the cn so- 
lutions are anti-symmetric, we expect that for even num- 
bers of peaks the chemical potentials of the two solution 
types should be very close. For odd numbers of peaks 
there are no cn solutions. 

Removing all factors in common to the two chemical 
potentials we find from Eq. (^) and Eq. (Fn for the cn 
solution and from Eq. (|2^) and Eq. ( pif ) for the dn solu- 
tion: 



Mcn( m cn) OC 



2m r 



1 



E[m ca ] - (1 



i)K[m c 



Mdn(m-dn) OC 



(29) 



(30) 



where m cn is found by solving Eq. (^) numerically and 
mdn is found by solving Eq. (E4j) numerically. This gives 
the exact percent splitting: 



Q split 



= 100- 



Adn( m dn) - Acn(w cn ) 
■ |Adn( m dn) + Acn( TO cn) I 



(31) 



which falls off exponentially as m cril mdn —> 1 _ - 

In linear quantum mechanics it is known that for the 
spatial wavefunctions of two particles the antisymmet- 
ric superposition is higher in energy than the symmetric 
superposition. The exact opposite is true in the NLSE. 
Numerical studies of Eq. ([H]) show that the symmetric, 
nodeless, dn solutions are higher in chemical potential 
than the antisymmetric cn solutions. But just as in linear 
quantum mechanics, the more the wavefunctions overlap 
the higher the splitting. 



C. Complex symmetry-breaking solutions 



treatment is identical to that of the repulsive 



cas' 



By the same arguments used there, all intrin- 
sically complex solutions to Eq. (Q) may be written as a 
sum over standard elliptic integrals by the use of appro- 
priate Cayley transformationstS. 
Writing the wavefunction as: 

f{x) = r{i)e^ (i) (32) 

one may divide the NLSE into real and imaginary parts: 



{S'f = -2{±S 3 + ^S 2 - (3S + 2a 2 - 



a 
S 



(33) 



(34) 



where a and j3 are undetermined constants of integration 
and S = r(x) 2 is the single particle density |/(ai)| 2 . 



From substituting Eq. @ into Eqs. and 
again taking real and imaginary parts, four boundary 
conditions are obtained, the important one being phase 
quantization: 



(1) - 0(0) = 27m 



(35) 



where n is an integer which we will call the phase quan- 
tum number. 

Thus for the complex solutions, Eq. ( |33| ) and Eq. ( |3^ ) 
replace the NLSE as the equations to solve, together with 
four boundary conditions, of which phase quantization is 
the most i mportan t, an d the n ormalization, Eq. 



In Sees. IIIB1 and IIIB2 we showed that the real 



symmetry-breaking solutions have a density proportional 
to cn 2 and dn 2 , respectively. The former of these vanishes 
at 2j points around the ring, while the latter is nodeless. 
We look for complex solutions which effectively interpo- 
late between these two real solution types. The motiva- 
tion for such a solution will become clear in Sec. [VB. 
Using our physical intuition, we are again able to by- 
pass the use of Cayley transformations, as we did in the 
repulsive case. 

We remind the reader of the Jacobian elliptic function 
identity: 



1 



cn (x | to) = — (dn [x \ to) — (1 — to)) 



(36) 



We thus generalize the real symmetry-breaking solu- 
tions, Eq. (Q) and Eq. (Q, as follows: 



r 2 (x) = A 2 (dn 2 (kx + 8 | m) - 7(1 - to)) 



(37) 



When 7 = the dn solution is recovered; when 7 = 1 
the cn solution is recovered. As in the repulsive case, 
we temporarily set 6 — 0. j is to be interpreted as the 
number of peaks in the density r(x) 2 . We will consider 
the case of general S, and thus degeneracy, later. 

Using the solution methods as outlined in the repul- 
sive caseEj, k is set to 2jK(m) in order to satisfy the 
boundary conditions. The solution then becomes: 

r 2 {Z) = A 2 (dn 2 (2 jK(m)x \ to) - 7 (1 - m)) (38) 

The chemical potential and the parameters 7, A 2 , and 
a arc then: 



fj, = -- + l2j 2 \ 2 E{m)K(m) - 4j 2 {2 - m)X 2 K(m) 2 



-j~ 2 \~ 2 + 8E(m)K(m) 
8(1- m)K(m) 2 



A 2 = 8j 2 X 2 K{mf 



(39) 



(40) 



(41) 
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a = -L[ A~ 2 ((-l + 8j 2 \ 2 E(m)K(m) 
v2 

-8j 2 (l - m)X 2 K(m) 2 ) * (1 + 8j 2 A 2 is:(TO) 2 
+64/A 4 ^(TO) 2 J ft:(TO) 2 - 16£(m) (j 2 \ 2 K(m) 
+4/A 4 X(m) 3 )))] 1 /2 (42) 



This leaves the constant of integration a in Eq. (34) 



the interpolation parameter 7, the prefactor to the den- 
sity A 2 , and the chemical potential fx, determined up to 
the number of peaks j, the scale A, and the parameter 
to. For a given A and j we then numerically integrate 
Eq. ( |34| ) and adjust m until the phase is quantized on 
the ring, i.e. until n is an integer. 

All parameters are monotonic in to. Furthermore m 
has an extremely limited range. Because the solution d3§| ) 
was designed to interpolate between the two real solution 
types, outside of the small splitting between the value of 
to for the real cn and dn solutions (see Sec. Ill A3) there 
are no complex solutions at all. We find numerically that 
outside of these splittings a becomes purely imaginary. 
This would invalidate the hypothesis that the phase is 
real. Thus to find the complex solutions one need only 
scan the values of m between m cn and TO^n, so that the 
algorithm is quite straightforward. 

If there are an odd number of peaks then there is no 
cn solution. In this case we just consider the odd-peaked 
cn solution for the box, and treat it as a limiting case 
of the complex solutions. By symmetry of the ring, n 
can be either positive or negative, so that each solution 
is two-fold degenerate, just as we found for the constant 
amplitude solutions. 

In Fig. H we show the lowest two solutions. We have 
plotted the amplitude above the phase to make appar- 
ent that the phase is nearly constant over the peaks and 
jumps sharply between them. When the troughs go to 
zero, the phase becomes a step function and the height 
of each step approaches n, which recovers the cn-type so- 
lutions. As a consistency check, we note that all of the 
chemical potentials and other parameters in the com- 
plex case approached the values found in the real case as 
7 — > + and 7 — * 1 _ , respectively. 

If S is generalized so that it is arbitrary, a similiar 
degeneracy to what was found in Eqs. (21) and (28) is 
obtained: 



(43) 



where the factor of will be explained in Sec. [V A 



D. Spectra 

We show the chemical potential spectra as a function 
of A~ 2 for the four types of stationary states on the ring: 
real with nodes, constant amplitude, real without nodes, 
and intrinsically complex. In Fig. tfl the two lowest real 



spectra are shown. For comparison we have overlayed 
the four lowest constant amplitude spectra on the same 
figure. In Fig. |^ we show the spectra of the nodeless 
solutions: the three lowest for the complex ones and the 
five lowest for the real ones. 

The spectra in Fig. |^ are very linear because they are in 
the bright soliton regime. As soon as the peaks overlap 
appreciably these two solution types no longer become 
available. However, for the cn spectra shown in in Fig. || 
both the bright soliton limit, to the far right, and the 
particle-on-a-ring limit, to the far left, hold. 

For the experimentally reasonable test scale of A = 
1/25 the ground state is the single-peaked, nodeless dn 
solution. The first excitation is the two-peaked anti- 
symmetric cn solution, the second is the two-peaked dn 
solution, the third is the three-peaked n — ±1 complex 
solution, the fourth is the three-peaked dn solution, and 
so forth. 

Since the real solutions are a limiting cases of the com- 
plex solutions, the three symmetry-breaking types scale 
in the same way and their energy levels do not cross, as 
will be further explained in the following section. But 
the constant amplitude solutions depend differently on 
scale, so their energy levels can cross with those of the 
other solutions. 



IV. BOUNDS AND SYMMETRIES ON THE RING 

The stationary states presented above are the 
bounded, quantized version of bright soliton-trains. So- 
lutions on the ring have many special properties and sym- 
metries. We detail their bounds, degeneracies, and point 
group symmetry in light of the BEC. That is, we con- 
sider the effect of changing the scale parameter A, which 
corresponds to accretion of atoms into a condensate. Ex- 
perimental predictions will be made in Sec. 



A. Bounds 

Box-type stationary states can have an arbitrarily large 
number of nodes. But the excitation level of nodeless so- 
lutions is limited. We set three bounds on the nodeless 
stationary states: the maximum chemical potential; the 
minimum and maximum phase quantum number; and 
the minimum scale to obtain j notches. As a conse- 
quence of these bounds there are some scales at which 
no nodeless solutions exist. 

Mathematically, the maximum number of peaks that 
can fit on the ring is obtained from the lower limit on 
the period of the dn function in the solution (|||). When 
the peaks overlap too much they are no longer solutions 
to the NLSE. The dn function approaches its minimum 
period of n as to — > + . In this limit Eq. (|39|) becomes: 



jj-max — ~ ^ 2j 2 A 2 7T 2 



(44) 
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In this same limit the amplitude approaches a con- 
stant which the normalization constrains to be 1. From 
Eqs. (|34|) and ( f42| ) we find a relation between the maxi- 
mum number of peaks, the phase quantum number, and 
the scale A: 



A 



2 - 8n 2 

max ^ 



(45) 



If 7i = the dn solution is recovered. Substituting 
Eq. into Eq. (Q and setting n = we find that 
P'max = — 1- Therefore the chemical potential of the 
plane wave solution is an upper bound for all nodeless 
solutions. 

We can set a bound on the phase quanta simply from 
the parameter range mentioned in Sec. [II C. Because the 
complex solutions interpolate between the cn and dn so- 
lutions: 



< n < 



J 



(46) 



Note that n is stricly less than j/2. This is because the 
phase quanta of the cn solution is actually 0, not j/2. 
j/2 is the limiting case, and j can be odd, even though 
an odd-peaked cn solution does not solve the NLSE on 
the ring. 

In Fig. U we plot the inverse scale at which each j 
becomes available for the maximum n within the range 
of n. For odd j this is only half a phase quantum away 
from j/2, while for even j it is a whole phase quantum 
away. Thus we expect that the odd j solutions should 
become available at lower inverse scales than the even j 
solutions, due to the quantization condition: 



(A ) m 
(A" 1 ), 



tV-8 + : 



ttV-2 + Wo, 



(47) 



for even or odd j, respectively. The order of j-peaked 
solutions turning on as a function of A -1 is then: 
1,2,3,5,4,7,9,6,11,13,8,... Thus for large A" 1 many 
more odd than even solutions are available, as again may 
be seen in the plot. 

In Fig. |^ and in Eq. ( |47| ) it is apparent that the min- 
imum inverse scale for a complex solution is ny/2. This 
means that the ground state for less than about 4.5 heal- 
ing lengths to the box length is the n — constant so- 
lution, rather than the single-peaked n = dn solution. 
This is also the source of the constant factor in the en- 
tropy of the nodeless solutions, Eqs. (|2|) and (p3|). ir^/2 
is the minimum size of a peak. 

Supposing that A -1 is large there are then three 
regimes: for small j all solution types are available; for 
intermediate j even-peaked solutions are cn and go to 
zero, while odd-peaked solutions are complex and shal- 
low; and for large j only cn solutions are available, so 
that there are no odd-peaked solutions. 



B. Theory of point groups 

From Eq. ([j]) one may consider the negative nonlinear 
term in the NLSE as an effective potential generated by 
the wavefunction. As long as the peaks in the wavefunc- 
tion are well separated then the self-generated troughs do 
not interact. On a ring one may treat the well-separated 
limit as a rotationally symmetric set of j potential wells, 
where j is, as before, the number of peaks in the magni- 
tude of the wavefunction. In Fig. ^ the case for j = 5 is 
shown. 

This is in complete analogy to the quantum dynamics 
of a particle in a j-fold rotationally symmetric potential 
in two dimensions. Such a physical situation is described 
by the Cj symmetry point group. The number of irre- 
ducible representations of the Cj group determines the 
degeneracy of the solutions. Thus one may look up the 
Cj character tablesEa in order to find the degeneracy. 

Consider the case of j = 5: from Eq. ( p7| ) it may be 
seen that there is a real, n — 0, dn solution and two 
complex solutions, n = 1 and n — 2. As the two com- 
plex solutions are each two-fold degenerate there are five 
solutionsJji all. This is identical to the the C5 charac- 
ter tableEj: there are one real solution and two doubly- 
degenerate, intrinsically complex solutions. The latter 
are essentially a linear combination of independent bright 
solitons with coefficients which are simply the appropri- 
ate group characters, as familiar to chemists in molecular 
orbital theoryEJ. 

In the limit that A — > + , i.e. that the peaks are well 
separated, given j peaks there are j nearly degenerate 
solutions. If j is even there are 2 real solutions, one 
symmetric and one antisymmetric, and (j — 2)/2 complex 
solutionSjeach two-fold degenerate. The splitting is given 
by Eq. (|3l|). If j is odd there is only one real solution, 
the symmetric one, and (j — 1)/2 complex solutions, each 
two-fold degenerate. 

An upper bound may be placed on the splitting by ap- 
plying Eq. (|3l"| ) and keeping in mind that the odd-peaked 
cn solution is a limiting case on the ring. For both even 
and odd j the ordering of the chemical potentials is as de- 
tailed in Sec. |||d, with the anti-symmetric solution be- 
ing the lowest and the other symmetry-breaking solutions 
being higher in chemical potential as n — > + . There is 
an additional degeneracy due to symmetry-breaking, as 
seen in Eqs. @, (||), and @. 

Thus we have shown by group theoretic considera- 
tions that our bright stationary state solutions to the 
NLSE under periodic boundary conditions include all 
symmetry-breaking eigenstates of evenly spaced peaks 
identical up to a phase. 



V. EXPERIMENTAL PREDICTIONS 
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A. Formation of the ground state 

The symmetry of the ground state of the NLSE on 
a ring is scale-dependent. For A -1 < \2it the ground 
state is a constant; for A -1 > \/2it it is single-peaked, 
symmetry-breaking, and nodeless. The width of this 
peak is about \/2nX. Recalling the definition of A, for 
fixed L, i.e. a fixed trap size, we may consider the 
symmetry-breaking as a function of either particle num- 
ber or scattering length. The former case corresponds to 
the experimental situation of adiabatic accretion of par- 
ticles in quasi-lD. The latter corresponds to the tuning 
of a scattering length, by a Feshbach resonanceE^l or other 
means. 

Thus one expects under adiabatic growth of the con- 
densate to observe the density of the attractive BEC, as a 
function of the number of atoms trapped, to be constant, 
then to peak in the middle, thereby breaking symmetry, 
and finally to disappear from view as the width of the 
peak becomes less than the wavelength of the imaging ra- 
diation. In Fig. m we show this sequence of events in four 
stages superimposed on each other. Under box boundary 
conditions one expects a similiar occurrence. Although 
the ground state is not a constant, it is quite broad for 
A -1 ~ 1. As A -1 increases, the ground state becomes 
steadily sharper until it is no longer visible under the 
imaging radiation. 

Though it has been shown that the NLSE models the 
BEC quite well in the repulsive case, it has been unclear 
as to how well it models the physics of the attractive 
case. There have been relatively few experiments with 
attractive interactions, and analysis of the NLSE in a 
three-dimensional trapped condensate suggested a small 
upper limit on the number of. atoms that can be main- 
tained in a stable condensate^ Ho. 

We predict that in quasi-lD, i.e. for transverse dimen- 
sions on the order of £, a slow enough accretion of parti- 
cles to the attractive condensate will have the sequence 
shown in Fig. || for the ring and a similiar sequence for 
the box. An experiment looking for such a change in the 
shape of the ground state could answer two questions at 
once: whether there are bright solitons and whether the 
NLSE is a useful model for the attractive BEC. If it can 
be shown that the NLSE models the attractive BEC in 
quasi- ID then there are a plethora of rich phenomena in 
fiber optics which could have a direct analogue in the 
BEC. 



B. Modulational instability 

For A" 1 < V2tt we have shown that the constant solu- 
tion is the ground state. We have discussed the case of 
adiabatic changes in the particle number or the scattering 
length. Let us now consider non- adiabatic changes. 

We first note that the depth of the dn solutions is 
highly scale dependent. We solve Eq. ( ^4|) in the limit as 



to — > + for A 1 : 

A" 1 = V27rj (48) 

For each j G {1, 2, 3, ...} there is a scale for which the real, 
j-peaked dn solution is in fact a constant. Near these 
scales the dn has shallow modulations. Midway between 
these scales it is exponentially close to zero between the 
peaks, in Fig. ^ we plot such a sequence for the j = 5 
solution. The four plots of the wavefunction maintain the 
same phase quantum number, n = 0, yet vary greatly 
in depth with only a small fluctuation of scale. This 
extreme sensitivity in the form of the solutions to an 
external parameter, here A -1 , is akin to the modulational 
instability found in fiber optics, in which a continuous 
wave laser picks up highly variable sidebands in response 
to small changes in the powerQ. 

If the scattering length is positive then the BEC is re- 
pulsive and the ground state on a ring is the constant so- 
lution. Using a Feshbach resonance the scattering length 
may be tuned rapidly negative. In the case of 23 Na the 
condensate has been experimentally observed to disap- 
pear pfcam view when the scattering length becomes neg- 
ativeE3E3, with various ex&Laaations based on atomic re- 
combination mechanismsS-TEa. In the case of 85 Rb the 
condensate shrank to below the resolution limit and then 
emitted a burst of high-energy atomsE3. 

We have shown here that the stationary NLSE pro- 
vides a mechanism for the apparent destruction of the 
condensate even in the absence of recombination, at least 
in the quasi-lD approximation. For a fast transition in 
the scattering length we expect the constant solution, in 
response to noise, to be able to easily shift to a many- 
peaked soliton solution without having to change phase 
quantum number. It may even be possible to control the 
form of this transition by using intentional noise to in- 
duce the desired j-peaked soliton train solution. Induced 
modulation has proven quite useful in fiber optics, as for 
example in the creation of a new kind of laseitl. 

VI. CONCLUSION 

We have presented the complete set of stationary so- 
lutions to the nonlinear Schrodinger equation under pe- 
riodic and box boundary conditions in one dimension for 
the case of attractive nonlinearity. In a box all solutions 
may be taken to be real. On a ring there are four solu- 
tion types: constant amplitude solutions which are plane 
waves; real symmetry-breaking solutions with nodes; real 
nodeless symmetry-breaking solutions; and a novel class 
of complex, symmetry-breaking, nodeless solutions. All 
symmetry-breaking stationary states on the ring are de- 
scribed by the theory of point symmetry groups. 

Constant amplitude solutions and real solutions with 
nodes are in one-to-one correspondence with those of the 
analogous particle-on-a-ring and particle-in-a-box prob- 
lems in linear quantum mechanics. Nodeless, symmetry- 
breaking solutions are uniquely nonlinear. Solutions of 
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non-constant amplitude may be treated as bright soliton- 
trains. As the natural size of a density-peak is the 
minimum scale size needed to obtain nodeless solutions 
is X/£ = 7T-\/2. As a consequence, the form of the ground 
state depends on L/£, which in the context of the BEC 
means that as atoms accrete the ground state on the ring 
breaks symmetry. 

In addition to describing the properties and physical 
meaning of stationary states in detail, we have made ex- 
perimental predictions specific to the BEC. The creation 
of an attractive BEC in the quasi-one-dimensional regime 
could answer two questions at once: is the NLSE with 
attractive nonlinearity a good model; and can one ob- 
serve solitons? Thus far experiments have only been per- 
formed in the three-dimensional regime. The quasi-one- 
dimensional solutions presented in this work may suggest 
further experiments. -Elsewhere we have illustrated such 
solutions numerically!!!. 
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FIG. 1. This graphical solution of Eq. ([?]) shows that for 
a given scale and number of nodes the real solution to the 
stationary NLSE under box or periodic boundary conditions 
is unique. A is the scale and j — 1 with j 6 {1, 2, 3, ...} or j 
with j € {2, 4, 6, ...} is the number of nodes, respectively. The 
three curved lines are plots of Eq. ((?]) solved for the number 
of nodes j, with A -1 = L/£ = 10,25,50. The left-hand side 
of the plot is the m = 0, linear limit, while the right-hand 
side exponentially approaches the m = 1, bright soliton limit. 
The solutions are found where these lines intersect with the 
horizontal lines of j. Note the rapid convergence to m — in 
the high j limit, so that for large j the solutions are in the 
linear regime. 
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FIG. 2. Real stationary solutions to the NLSE under box 
and periodic boundary conditions. The solid lines show so- 
lutions in one-to-one correspondence with those of the analo- 
gous particle-in-a-box and particle-on-a-ring problems in lin- 
ear quantum mechanics. The dashed lines show uniquely non- 
linear, nodeless solutions found only on the ring. Both solu- 
tion-types may be characterized as bright soliton trains. Box: 
the solid lines in (a)-(d) are the ground state and first three 
excited states. Ring: (a) nodeless, symmetry-breaking ground 
state (note that the solid line has overwritten the dashed one) 
(b)-(d) symmetric and anti-symmetric solutions (c) only the 
dashed line is a solution. All plots are for the test scale of 
£/L = 1/25. 



FIG. 4. Two intrinsically complex, stationary, bright soli- 
ton train solutions on a ring, j is the number of peaks, n 
is the phase quantum number, and all plots are for the test 
scale of £/L — 1/25. (a) Amplitude and (b) phase/27r of the 
j = 3, n = 1 solution, (c) Amplitude and (d) phase/27r of the 
j = 4, n = 1 solution. Note that these solutions are two-fold 
degenerate, as the chemical potential depends on n 2 . 
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FIG. 3. Chemical potential spectra of real stationary 
states with nodes, as a function of inverse scale L/£, with 
stationary plane wave spectra shown for comparison. Dashed 
lines: shown are n = 0, 1, 2, 3, where n is the phase quan- 
tum-number of the plane wave on the ring. Solid lines: real 
stationary states of the NLSE in a box and on a ring are soli- 
ton-trains. Shown are j = 1, 2, 3, 4 with j — 1 the number of 
nodes in a box and j = 2, 4 the number of nodes on a ring. The 
linear regime to the far right corresponds to the bright soli- 
ton limit in which the peaks are well separated, while the far 
left corresponds to the particle-in-a-box or particle-on-a-ring 
limit. 
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FIG. 5. Chemical potential spectra for symme- 

try-breaking, nodeless solutions on the ring. These solutions 
have no analogue with the particle-on-a-ring solutions in lin- 
ear quantum mechanics, j is the number of peaks, n is the 
phase quantum number and L/£ is the number of healing 
lengths per box length. Dashed lines: j = 1,2,3,4,5 from 
left to right. Solid lines: (j,n) = (3, 1), (4, 1), (5, 1) from 
left to right. The spectra are nearly linear because when the 
peaks overlap appreciably they no longer solve the nonlinear 
Schrodinger equation, so that they only exist in the bright 
soliton regime. 
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FIG. 6. Minimum inverse scale for j peaks to become 
available. The lower curve is odd j; the upper curve is even 
j. Note that in general there are many more odd solutions 
than even solutions available. The ordering of the solutions 
is j = (1,2, 3, 5, 4, 7,9,6,11,13,...). 




FIG. 7. Five-fold group symmetry on the ring. The dips 
are the mean-field effective potential produced by the five 
peaks in the condensate. 
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FIG. 9. Modulational instability: fluctuations in par- 
ticle number lead to symmetry-breaking solutions of radi- 
cally varying depth. All four plots are for j — 5 peaked 
dn solutions, with scale varied from the minimal num- 
ber of healing lengths needed to get a five-peaked solu- 
tion to just over the minimal number needed to obtain a 
six-peaked solution, (a) L/£ = 5tv\/2, jl — —1.0000 (b) 
L/i = 5tt\/2 + 10~ 2 , jl = -1.0028 (c) L/£ = 25, jl = -1.4657 
(d) L/i = 671-^/2 + 10~ 3 , ji = -1.0001 
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FIG. 8. Adiabiatic formation of the ground state 
on a ring: the onset of symmetry-breaking. From 
the broadest to the thinnest solution the scales are 
L/£ = 7r\/2, 5, 10, 25. The respective chemical potentials are 
jl = -1, -1.852, -6.250, -39.025. In a box the formation is 
similiar, except that at the broadest scales the ground state 
goes smoothly to zero rather than being constant. 
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